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ABSTRACT 

We extend QCD sum rule analysis to moderate energy fixed an- 
gle Compton scattering. In this kinematic region there is a strong 
similarity to the sum rule treatment of electromagnetic form factors, 
although the four-point amplitude requires a modification of the Borel 
transform. To illustrate our method, we derive the sum rules for helic- 
ity amplitudes in pion Compton scattering and estimate their large-t 
behavior in the local duality approximation. 
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1 Introduction 



The method of QCD sum rules originated from an analysis of 
the high energy behavior of two-point correlation functions of QCD 
currents. More than 10 years ago, this method was extended to elastic 
form factors |2|, [3| by treating three-point functions. In the sum rule 
method, one considers Green functions of interpolating fields given 
by current-correlators which have nonzero resonant contribution on 
selected lowest states. The non-perturbative resonant region is then 
enhanced over the continuum by using Borel transforms p| |3[. A 
problem in the more general cases M involving four or higher-point 
functions is the larger number of invariants on which the function 
depends. In this paper, our goal is to extend the QCD sum rule 
approach to fixed angle pion Compton scattering, a process described 
by an amplitude which is closely related to the pion electromagnetic 
form factor. 

Our considerations build on previous work on sum rules for the 
pion form factor in Refs. || and ||, and are made possible by the 
observation that, at high energy and fixed angle, the underlying quark- 
photon scattering subprocess has a short-distance nature. 

The approach that we investigate is also related to the observation 
(see, e.g., f|) that QCD sum rules can be formulated as a specific 
variant of the finite-energy sum rules analogous to those applied orig- 
inally to high-energy hadron-hadron scattering || 0, |§. Compared 
to the sum rule treatment of form factors 0, we must take into 
account the presence of additional cuts in the w-channel when we con- 
struct a dispersion relation. These additional cuts force us to work 
in a kinematical regime in s and t where the cuts are far from the 
region enclosed by the integration contour. Our aim is to show that 
if we choose s, t and s + t to be moderately large in the physical 
region s > 0,t < 0, then a dispersion relation can be written down 
and, hence, a QCD sum rule can be derived. These sum rules re- 
sult in invariant amplitudes with s and t dependence in general. This 
dependence is not expected to be the same as for perturbative calcu- 
lations based on elastic scattering factorization theorems 0, because 
the sum rule calculation begins at 0{a° s ). Thus, when our method 
is extended to the proton, it may be possible to confront both the 
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sum rule and the standard perturbative calculations with experiment 
at a fixed energy. This would help clarify some controversial issues 
concerning the relative importance of these contributions at available 
energies 



10, 11, 12 



Our arguments are illustrated by an explicit calculation of the 
discontinuity in the perturbative amplitude at the lowest order, and 
by a derivation of the sum rules for invariant amplitudes in the case 
of Compton scattering of the pion. 

In Section 2, we reformulate the QCD sum rule approach to the 
pion electromagnetic form factor. To this end, we modify the standard 
procedure of Refs. [I], |], [| by employing an integral representation of 
the Borel transform. In Section 3, we develop the relevant dispersion 
relation for the four-point amplitudes, and derive a sum rule by using 
a variation of the Borel transform. In Section 4, we discuss the iso- 
lation of the pion contribution. Section 5 contains our results for the 
lowest order spectral functions and the QCD sum rule estimate for 
the Compton amplitude in the "local duality" approximation. Our 
conclusions are summarized in Section 6. 

2 Form Factor Sum Rules 

The most popular choice of the interpolating fields in the pion case is 
the axial current 

= u(x)7^ 5 d(x), (1) 
whose projection onto a single-pion state \p) is given by 

(®K{x)\v)=ihv»e-^ x . (2) 

The electromagnetic current J M defining the pion form factor is given 
by 

2 1 - 

Jn = ^u^u - -d^d. (3) 

To study the form factor, one analyzes the 3-point correlation func- 
tion of the electromagnetic and two axial currents |2], [3| 

rW(pi,p 2 ) = f(0\T ( Va (y)J u (0) V l(x)) |0)e-^ + ^d 4 :rdV (4) 
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Analogously, to investigate pion Compton scattering, one should 
consider the 4-point function involving two axial and two electromag- 
netic currents 

I\,a^(Pi,P2,M) = i j ' d 4 xd 4 yd 4 ze- l ^ x+ip2 -y~ l ^ z 

x{0\T(pMJ v (z)J x (0)r,l(x))\Q), (5) 

where the photon lines are assumed to be on-shell. 

Let T 3 (p1,p2,t), with t = (p 2 — Pi) 2 , denote one of the scalar 
invariant amplitudes isolated from the tensor expansion of T^ u \ |fj. |3| . 
The amplitude T 3 obeys the double dispersion relation 

T 3 (p?,»2,t)= / ds% [ ds 2 -, 3 ^,1' 2 ' - , + subtractions, (6) 

Jo Jo (si -Pi)(s 2 -pi) 

where A 3 (si, s 2 , t) is the double discontinuity of T 3 across the branch 
cut associated with thresholds in p\ and p\. 

The next step within the QCD sum rule approach is to introduce 
the double "Borel" transform $ 3 (Af?, Mf, t) |, | of T 3 (pf,p|, t). The 
Borel transformation can be defined in different ways. For our pur- 
poses, the most convenient is the integral representation 

*3«M|, ( ) = _^/ c «/ o | e - m « r3(p ^, () , (7) 

where C denotes any contour that encloses the branch cuts in the p\ 
and p 2 planes, from zero to infinity. For any non-zero values of the 
Mf, the contours C may be thought of as closed across the cut at 
positive infinity, where the integrand is exponentially suppressed. 

For the 3-point function, the integral form fl7|) is completely equiv- 
alent to the standard procedure P], |^, [| for constructing the Borel 
transform. However, as we will see, it is very suggestive of an exten- 
sion to the four-point function related to fixed-angle high-energy pion 
Compton scattering. 

Note that eq. (|7|) immediately gives zero result for the subtraction 
terms in T 3 {p\,p\,t), since they are polynomials in the pf. Applying 
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(0) to (§ then gives 



* 3 (Mf , M|, t) = f ^ f f| t)e -./-f-«. (8) 

To derive a sum rule, we employ our freedom in choosing the origi- 
nal contours C in the definition of the double transform, eq. (0) . Two 
relevant choices are shown in Fig. 1. Contour C a encloses the branch 
cut closely, while Cb encloses it with a circle of radius ( 2 about the 
origin, and joins with C a for \p 2 \ > ( 2 . Let us suppose that we choose 

C 2 » ^ (9) 

where /i is the scale associated with nonperturbative power correc- 
tions in T^,{p\,p\,t). The inequality makes it reasonable to apply the 
operator product expansion along Cb- Around the contour Cb, we may 
now approximate 

T^{p\,p\,i) = T pert {p\,p\,t) + power corrections, (10) 

where T^ pert {p\, p\, t) is the perturbative expansion for T3 starting with 
the diagrams of Fig. 2. The power corrections indicated in eq. ( |10"D 
are proportional to the familiar vacuum condensates, 

(Ol ^g^ggp) {0\a s (qg ) 2 \0) 

(Pi) 2 2 (p?) 

with Ki and K 2 being, in general, functions of t/pf and p\jp\ ■ We 
shall assume that ( 2 = 0{t) in eq. (§), so that possible higher-order 
logarithmic corrections in these variables are not large. 

Around the contour C a , the approximation of eq. (|10"D is clearly 
inappropriate. Instead, for small values of p\ and p\, we parame- 
terize T 3 phenomenologically - in terms of dominant resonances plus 
continuum: 

A 3 (p 2 1 ,p 2 2 ,t) = Af en (p 2 1 ,p 2 2 ,t). (12) 

The sum rules now result from eq. (H) by equating the two transforms, 
one around C a and one around Cb, using the two approximations, eqs. 
and (|12"D respectively, for A 3 , 

_ r°° d,Si f°° ds 2 , p er t /,phens - Sl /Ml~s 2 /Mi 

° " Jo WJo M? (As ~ As )C 



power corrections = A 1 h A 2 -r-^ ( 1 1 ) 



+ power corrections, (13) 
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with A^ pert the double discontinuity of T^ pert around the p\ and the 
p\ cuts. Here the contributions of power corrections will appear only 
from the (counterclockwise) circle of radius ( 2 for C a in Fig. 1, 

J_ f dp 2 e-v 2 l M2 = (-1)" 

27riV|=e 2 M 2 {p 2 ) n " (n- 1)!(M 2 )™' 1 j 

The remainder of the C a integral gives zero for power corrections, since 
they are well approximated by functions that are analytic except at 
the origin. Of course, this is an approximation, but an adequate one 
for ( 2 ^> /i 2 , as specified in eq. (fj|). The power corrections have been 
computed in Refs. 0, [|. Once the power corrections are specified in 
eq. fllTf ), information can be extracted when an explicit form is given 
for A^ phen . For the form factor, the standard choice "first resonance 
plus continuum" is made in Refs. @, U : 

A phen ( Sl ,s 2 ,t) = T^t)6( Sl -m n 2 )6(s 2 -m n 2 ) 

+ Ar rt (s u s 2 ,t)(l - 6(s - Sl )9(s - s 2 )), (15) 

where T£(t) is the pionic scalar amplitude and s is a parameter spec- 
ifying, in this model, the onset of continuum. Using (|T5| ) in (|T3|), we 
find the M 2 - and s - dependent sum rule, 

Tg(t) = [ S °d Sl [ S °ds 2 Ar t (s 1 ,s 2 ,t)e-^ +s ^ M2 
Jo Jo 

+ power corrections, (16) 

where we have set M 2 = M\ = M 2 and neglected m 2 compared to all 
other scales. 

Note that contributions from the states beyond y/so in mass have 
cancelled here. Fixing so from the requirement of the best stability 
of the right hand side of eq. (0) with respect to variations of the 
Borel parameter M 2 (this gives so ~ 0.7 GeV 2 ) and then taking the 
limit M 2 — > 00 (which eliminates power corrections), one arrives at 
the "local duality" formula for the pion form factor || . We now turn 
to the application of this reasoning to the four-point amplitude. 
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3 Sum Rules for Four-Point Amplitudes 



In general, the analytic structure of a four particle scattering ampli- 
tude, or a four-current correlation function, is much more complicated 
than that of a three-point function. Nevertheless, at high energy and 
momentum transfer, the space-time structure of the Compton scat- 
tering amplitudes simplifies somewhat, and becomes similar to the 
corresponding form factor amplitudes. 

Consider, for instance, the four-current Green function T u \^ a , eq. 
@. The aim of the sum rule approach is to derive information about 
F u \fia at vl — m 1i where this function develops bound-state poles. Of 
course, we cannot compute Y perturbatively at these momenta. If, 
however, we can use a dispersion representation to relate the Green 
function in this nonperturbative region to its behavior in a region 
where perturbation theory and the operator product expansion may 
be applied, the road is open to apply the methods described in the 
previous section. We shall argue that it is possible to do this quite 
directly when the Mandelstam invariants s = {pi + qi) 2 and t = 
(Pi ~ P2) 2 in eq. (|5|) are large, with the ratio s/t fixed. 

The first step in this process is to note that even for on-shell ex- 
ternal particles, the underlying hard scattering at fixed s/t remains 
a short- distance process in any two-to-two scattering amplitude [|T3 



Perturbation theory runs into trouble from virtual particles that ei- 
ther (1) are on-shell and parallel to one of the four external momenta 
(collinear lines), or (2) have vanishing momenta in all four compo- 
nents (soft lines). Momentum configurations of this type will typi- 
cally diverge in the zero- mass limit, or, equivalently, suffer from large 
logarithms of ratios of external energies to masses. The situation 
for Compton scattering at lowest order in electromagnetism, however, 
simplifies when the external photons carry physical polarizations. In 
this case, virtual lines parallel to the incoming or outgoing photons 
do not produce large logarithms, nor do soft lines attached to such 
collinear lines. The failure of perturbation theory, then, results from 
the lines that are collinear to one of the hadron momenta and/or 
from the soft quanta that connect them. This is qualitatively the 
same situation as for the elastic form factor. Indeed, in space-time, 
the incoming photon and the outgoing photon originate at the same 
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point. This observation has been used elsewhere in the context of a 
perturbative QCD calculation of Compton scattering to show that the 
phases of this amplitude are free of infrared divergences |TJ . 



To take advantage of these observations, we shall need to develop 
the relevant dispersion relation. Its form will be slightly different 
from the form factor case, requiring the modified version of the Borel 
transform mentioned above. Let us describe how this works. 

Let T^(pf,pl, s, t) denote a particular scalar amplitude derived 
from T u \^ a in eq. fl5|). We will give a specific example below. Mo- 
mentum assignments are illustrated in Fig. 3. We will be interested 
in s and t large compared to m 2 , with —t/s fixed but not equal to 1 
or 0. We take qf = gf = 0, appropriate to physical photons. 

In general, we are not able to write for T 4 a double dispersion 
integral of precisely the form of eq. (0). This is because, when p\ or 
p\ becomes large enough, the kinematic identity 

s + t + u = pl+p 2 2 (17) 

forces u to threshold for fixed s and t. This means, first, that branch 
points in u appear at p\— dependent positions in the p\ plane (and 
vice- versa). Second, it means that, even at large p\-\-p\, a phenomeno- 
logical approximation like ([To]) would need to include resonances in the 
li-channel. Fortunately, this is not the end of the story. It is still pos- 
sible to write a useful double dispersion relation for T4. This relation 
employs the contour 7 shown in Fig. 4. 

Instead of running to infinity as in eq. 7 is closed in a circle 
of radius A 2 , where A 2 is a constant which we may take to satisfy 

A 2 = (s + t)/4. (18) 

With this choice, the w-channel threshold never comes closer than 
(s + t)/2 to 7 in the p 2 plane, when p\ is chosen to run along the 
corresponding contour in its complex plane. We therefore have, in 
place of eq. (|6|), 

rp / 2 2 ±\ 1 f . f . T i ( Sl ,S 2 ,S,t) 

We note that the portion of 7 that runs along the real axis gives an 
integral over a discontinuity, as in (0), but, generally speaking, the 
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discontinuity is not purely imaginary. Hence, the amplitude we are 
interested in will be complex in general, even for moderate negative 
p\ and p\. It remains, however, an analytic function of s if we are far 
from its thresholds. 

More formal arguments for the validity of this dispersion relation, 
and for the applicability of perturbation theory to the amplitude at 
large values of \pf\ are given in Appendix A. Here we shall assume 
that these results hold, and show how they lead to the desired sum 
rule. 

To derive a sum rule, we would like to apply the same Borel trans- 
form, eq. (0) to T 4 in eq. (|i"9|) that we applied to T 3 in eq. (H). The 
contours C in eq. (|7|), however, enclose the entire real axis, and would 
hence cross the contour 7 of eq. ([19]). This makes eq. (|7|) unsuitable as 
it stands. This technical difficulty can be circumvented if we choose 
a slightly modified version of the transform. As shown in Fig. 4, let 
T be any counterclockwise contour, enclosing the positive real axis of 
the pf plane, which starts at A 2 — ie and ends at A 2 + ie, but remains 
inside the overall contour 7 of Fig. 4. Note that T is not quite closed 
on the real axis. We assume that s and t are large enough so that we 
may choose 

A 2 > M 2 > /i 2 , (20) 

with the scale of power corrections, A 2 given by eq. (0) and Mj 
the arguments of the transform. We then define our modified "Borel" 
transform of T4 by 

41 1 ' 2 ' ' } ~ (2tt2) 2 Jr Ml Jr Ml ^ ) 

xe -*i/*S (l - e-^- p ^) T 4 (plpls,t) 

x (l - e -^- s ^ M2 ) (l - e -(^)/Mf ) . (21) 

In the second equality, the evaluation of the pj integrals has been car- 
ried out, using Cauchy's theorem and the representation eq. ([19]) for 
T4, by simply shrinking the V contours to the fixed points A 2 ±ie. For 
Si real and < A 2 , we pick up a pole contribution, just as in eq. ([8]) 
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for the three-point function, while for S{ on the circular part of the 7 
contour in Fig. 4, we get no contribution, as the T contour shrinks to 
the infinitesimal segment between ±ze without encountering any sin- 
gularities. At the special points, Sj = A 2 , the factor 1 — e-( x2 -"p\)l M 't 
vanishes at the pf = A 2 pole, so that the transform gets a vanishing 
contribution from these points as well. As usual, A 4 is the double 



discontinuity of T 4 . From eq. (j21|) , we will derive sum rules exactly 



analogous to the form factor result, eq. ([13]), by equating the per- 
turbative plus power corrections expression for T 4 (si, s 2 , s, t) with a 
phenomenological parameterization of A 4 in eq. (19). Let us see how 
this works for the case at hand, the Compton scattering of a pion. 

4 The Pion Compton Amplitude 

The pion matrix element 

M yX = ij d A ye-^(p 2 \T (J„(y)J x (0)) \ Pl ) (22) 

is the object we will study in terms of its invariant amplitudes, for 
which we may construct sum rules as in the previous section. M v \ 
can be expressed in terms of two invariant amplitudes. We use helicity 
form factors Hi and H 2 defined by the relation Ul5| 



M Am = Fie^V 1 ^ + F 2 e (2)A e^, (23) 

with 

tvtX pX 
e (l)A = ^ = e (2)A = (24) 



N" = ^P,r u R p P a =^ + ^-^ Pi - jR +7 2 - jR , (25) 

where 

v = Pl ' P2 - Sl (26) 

and where 

R = qi + q2 r = q 2 -qi. (27) 
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The "polarization" vectors e'*' satisfy 

e «. gi = e w -g 2 = 0, 
e« • e (i) = -<J y . (28) 

Note that these relations hold for all positive Si = p\ and s 2 = 
whether or not they are equal. We also note that e^ 1 ' serves as the 
normal to the scattering plane defined by the Compton process, while 
is a unit spacelike vector whose spatial projection is in the scat- 
tering plane, but which is orthogonal to the g^'s. Other expansions 
are also possible, but we choose this one for definiteness and for its 
simplicity. 

Now M v \ appears as the residue at the double pion pole at p\ = 
v\ — m 1 m the four-point Green function T u \^ a introduced in eq. 
To be more specific, let us define a double discontinuity A uXfia . of T u x^ a 
across the p\ and p\ cuts by 

Kx^ = r++ - r + _ - r_+ + r__ , (29) 

where we define T ++ = T u x f _ lCT (Pi + ^iV\ + it), and so on. The poles 
may be isolated from the x — > — oo, y — > oo limit of T v x^ a in eq. ([S]), 
by inserting complete sets of states in the resulting time ordering, as 

= jr^y J d 4 xd 4 yd 4 z d 4 p' \d 4 p' ' 2 9(y )9(-x Q ) 

x e — ipi-x—iqi-z+ip2-y 

x5(p'l-ml)5(p' 2 2 -ml)(0\ Va (y)\p> 2 ) 
x(^|T(J,(,)J A (0))b' 1 )(p' 1 |r / t( a; )|0) + I™* (30) 

where T cont includes contributions from the higher states and other 
time orderings, which do not contribute to the double pion pole. As 
usual, we express the matrix elements of the axial currents as 

(0MvM) = iUr^e-^^, (31) 

(pH( x )\°) = - i f*P'i/ P ' VX - (32) 
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This relates the single pion states p'i,p' 2 to the vacuum at times ±00, 
with f n the pion decay constant. We then derive the relation 

A^a m <t = fl V\pV2a (2tt) 2 5 {p\ - m\)b{p\ - m 2 ) 

x M vX (p 1 ,p 2 ,q 1 ) + A amt , (33) 

with M v \ given by eq. (|2~2"D , and with A cont the double discontinuity 

q£ pcorii 

We now relate the tensor A^o- to the invariant amplitudes Hi 
and H 2 of eq. fl23|). To do so, we begin by defining two projections, 

^ 12 \p 1 ,p 2 ,q 1 ) fl(T = -^T^, (34) 

r (1) (pi,p 2 ,giV = e w e w^. (35) 

When the pion poles in I^ 12 ) and T^ 1 ' are isolated, their coefficients 
will be proportional to Hi + H 2 and Hi, respectively. Thus, we want 
to identify a set of tensor structures in I^ 12 ) and I^ 1 ' in which the pion 
poles may appear. 

Note that A u \ flcr may be expanded in terms of tensors made from 
g a<3 and its external vectors, q, L and Pi . Following Refs. and @, 
we choose to expand in a basis made from the independent vectors 
Pi Pi — P2 an d qi- (The tensors ej^e^ built from the vectors in 

eq. (|2~3| ) may be expanded in this basis as well, when we include the 
metric tensor.) We then seek a vector n, which has the properties 

(n ■ qi) = (n ■ pi) — [n ■ p 2 ) = n 2 = 0. (36) 

It is easy to verify that no such n M exists with real coefficients. In fact 
we can immediately give an explicit expression for in terms of the 
polarization vectors above, as 

n»= [e^±ie^y . (37) 

The specific values of the n M are frame-dependent of course, and we 
need not give them here. We note, however, that either choice of 
sign in eq. fl3"TD will give the same (real) answers below. The overall 
normalization of also cancels in the sum rule result. 
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If we saturate the ji and a indices with an n that satisfies eq. ( |36| ) . 
we will project onto the tensor structure {j>\ + p 2 )n(pi + Vi)<j on ly: 
whose corresponding invariant functions we denote as 



Y { t 2 \ V ls,t) = nW( 12 )( Pl ,p 2 , gi ) MCT , (38) 
ri 1} (p?,M) = nW«( Pl ,p 2 , gi ) M(T . (39) 

Finally, taking as above the double discontinuities, denoted as A4 , 
we isolate the desired invariant functions 

A? 2) (p<,M) = f 2 An- Pl )(n-p 2 )(2n) 2 5(pl-ml)5(p 2 2 -ml) 
x (-g vX M vX ) + Af )cont 

= f%(n- Pi) (n • p 2 ) ( 27r ) 2 ^(P? - m l) s (pl ~ m l) 

x {H x {s t t) + H 2 (s,t)) +Ai 12)coni , (40) 

and similarly 

^\p 2 ,s,t) = fl{n.p 1 ){n.p 2 ){2^f5{p\-ml)5{pl-ml) 

x H^t) +A<i )cont . (41) 

The sum rules for the are derived as for T3 above. Consider, 
for instance, r^ 12 "*. First, we approximate by perturbation theory 
plus power corrections for a choice of contour 7 in eq. (0) that stays 
away from the real axis, except at p\ = A 2 , 

T^ 2 \pl, pi, s, t) = rji pert (Pi,pl, s, t) + power corrections. (42) 

This is analogous to eq. ([H]) for the three point function. Next, choos- 
ing 7 in eq. ( |i~9"D on the real axis, we approximate its discontinuity by 



a "first resonance plus continuum" form, analogous to eq. flT5|): 

Ai 12) = f 2 An.p 1 )(n.p 2 )(2n) 2 S(pj~ml)S(p 2 2 ~ml) 
x (H 1 (s,t)+H 2 (s,t)) 

+ Af )pert (1 - 6(s - Sl )e(s - s 2 )) . (43) 
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Combining these expressions in the transform of eq. fl2~T|), we derive 
the full sum rule as 



= r so d S1 r°ds 2 Af^( S1 , S2 , s ,t)e-^ 2 --/^ 

Jo Jo 

X (l- e -(A 2 - S i)/M^ ^ _ e -(X 2 -s 2 )/M^ 



+/l(s , t , ^ fry + /2(s , t , m?) <?!^!2> + .... (44 ) 



where we have written out the power corrections schematically, as 
in eq. (|TT|). The functions fi(s,t,M?) may be calculated as for the 
three-point function |2|, [|. We reserve this calculation for future in- 
vestigation. 

As in the form factor case of Section 2 above, So should be fixed 
by demanding best stability with respect to variations in M 2 = Mf = 
M|. We shall assume that the power corrections are small enough so 
that, as we increase the Mf, still keeping s/Mf large, there is a range 
where they are negligible compared to the perturbative contributions 
on the right-hand side of eq. (|44]). Then we may write a "local duality" 
sum rule, without power corrections: 

fl(n- Pl ){n-p 2 ) (H^t) + H 2 {s,t)) = f° da x P ^ 2 Af )peri ( Sl , s 2 , s, t) . 

Jo Jo 

(45) 

We should emphasize that the derivation of this result from eq. ([44]) 
involves assumptions on the behavior of the condensates. The right 
hand side of eq. (|^) is approximated only by the local duality integral, 
and is to be calculated using perturbation theory. This is the subject 
of the following section. 



5 Local Duality Sum Rules and Asymp- 
totic Behavior 

The lowest order perturbative spectral weights A^ 61 * can be com- 
puted from Figs. 5(a) and 5(b) by use of Cutkosky rules. The remain- 
ing lowest order diagrams, in which q 2 and p 2 are exchanged, do not 
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contribute to the double discontinuity in this kinematic region. The 
double discontinuity in p\ and pi is found by simply replacing the 
propagators that carry momenta (pi — k)^ and (p2 — by factors 
(2ii)5 + (k 2 ), and so on, where the plus indicates that energy flows in 
the same direction as for pi and p 2 . The procedure is similar to the 



one adopted in the case of the triangle diagram |2], [L6| for 3-point 
functions. 

We then write [\^ vert as the sum of the two spectral functions 
A$ pert and A^ peri found from the diagrams of Fig. 5, with an overall 
constant separated out, 

Ap ri ( Sl)S2 ,M) 

5 (A»^( Sl , S2 , S ,t) + Ap ri ( Sl , S2 , S ,t)). (46) 



3(2tt 

The extra factors are from the flavor structure of the current and 
the color traces, as well as from the factors of (27r) _1 associated with 
the loop integral. A factor of (27r) -2 associated with the dispersion 



relation eq. (|T^) has also been absorbed into the definition of the 
spectral function. 

For specific calculation, we found it useful to work in a "brick-wall" 
frame for pi andp2, withpi^ = P21. = 0; the corresponding kinematics 
is described in some detail in Appendix B. Let us consider the — g uX 

('12') 

projection, A4 , related by the sum rule to the combination Hi + H 2 . 
We have, in the normalization of eq. ([461), 



A (i2W = 16 1 d A k 



6 + (k 2 )6 + ((Pi -kf)5 + {{p 2 -kf) 



(Pi~k + qi) 2 
x (n • k) \(qi ■ pi + qi ■ p 2 ) (n ■ pi - n ■ k) 
— (qx ■ k)(n -pi +n -p 2 — 2n • k)], (47) 

with a related form for A| ) 12 ' )pert , found by replacing q\ by — q 2} and 
hence s by u in the resulting answers. For A^, we have 

A(l)pert = \\{l2)pert 
a 2 a 

+8 f d , h S + (k 2 )5 + (( Pl -k) 2 )5 + ((p 2 -k) 2 ) 
J (Pi~k + qi) 2 
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x (e^ ■ k) (n ■ k) |(e^ - A;) [n ■ pi + n ■ p 2 — 2k ■ n) 

+( e (i). n ) ((pi-Q-fa-k))}, (48) 

where we note that the second term is at least quadratic in k. The 
local duality sum rules are found by substituting the results of these 
integrals into eq. (|4"5]) . For instance, we find 



where 



f n 2 (H 1 + H 2 ) = dsj ds 2P T rt (s u s 2 ,s,t), (49) 
Jo Jo 



pTW,*) = 7 h ^ 2)pert ( Sl ,s,t). (50) 

{n-pi){n-p 2 ) 



The integrals in eqs. (|47|) and (^8|) are all elementary. The specific 
complete integrals and results are given in Appendix C. Here we shall 
only discuss the leading behavior at fixed angles. It is easy to show 
that it comes from the term linear inn - k, and is specified by 



10 (s x + s 2 \ ( (s-uf 



3(2tt) 2 V (~t) 2 J \ 2s(-u) J 

2p i l )vert . (51) 



(12)pert 

Pa 



Here we have used eq. ([50]), which defines p in terms of A, and have 
combined the contributions from both diagrams of Fig. 5. This re- 
sult exhibits a comparable power suppression to the spectral weight 
encountered in the form factor 0, [J . (Although the three point func- 
tion involves one less external photon, a factor of (p± +p 2 ) x is factored 
from its amplitude, so that the remaining form factor has the same 
scaling behavior as in Compton scattering.) To be specific, the lowest- 
order perturbative spectral weight in the latter case may be written 
as exactly 

Pvr3 ert (Sl,S 2 ,t) 

3 




t) 2 -4s 1 s 2 ) 1 /2- 



(52) 
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Asymptotically at large t this function has the characteristic behavior 



P^^(^^-Oil/h). (r )3) 



This is suppressed in t by one additional power compared to the pre- 
diction derived from the dimensional counting rules JT7], [18], [!|. 

The leading Compton scattering integral in eq. fl4~9|) can be evalu- 
ated trivially in the local duality approximation, in which we neglect 
Si/Mf compared to unity. We get, for H\ + H 2 , for instance, 

,„ T T \r2 10 f s ° , f so i (si + s 2 )A s 2 + u 2 



3(2tt) 2 Jo Jo (-t) 2 V 2su 



10 t S « \(x-t±±\. (54) 



3(2tt) 2 \(-t) 2 J\ 2su 

As in the case of the form factor, the result is suppressed by an ad- 
ditional power of t compared to the (1/i) perturbative behavior ex- 
pected from the dimensional counting rules. This may be understood 
in terms of the result that the true asymptotic behavior of form factors 
is dominated by hard gluon exchange ||. It is easy to see (Appendix 
C) that as s — > oo, with s\ and s 2 fixed, the momentum of the unscat- 
tered parton (k in Fig. 5) is forced to zero, while the incoming and 
outgoing momenta of the scattered parton approach the momenta of 
the external "pions" . This description of elastic scattering, mediated 
by a single parton that carries essentially all the energy of a hadron, is 



often called the "Feynman" mechanism ||19|| . Although it cannot give 
the true asymptotic behavior, it may be important at intermediate 
energies [10 



6 Conclusions 

In this paper, we have shown that QCD sum rules can be formulated 
for pion Compton scattering, generalizing previous results on meson 
form factors. We have shown that it is possible to write down a dis- 
persion representation of the scattering amplitude on a finite contour 
in the complex plane of the external masses, and we have explicitly 
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calculated the leading behavior of the spectral function. To derive 
a sum rule, we found it useful to introduce a modified definition of 
the Borel transform through an integral representation, eq. The 
sum rule, eq. (f44|), in the local duality limit, shows suppression by an 
additional power of t compared to the perturbative result, as shown 
in eq. fl5|). 

In the future, we hope to pursue the QCD sum rule procedure 
futher, with the inclusion of power corrections assoicated with the con- 
denstates that appear in power corrections from the operator product 
expansion. Here, however, we have demonstrated that such a program 
is possible. 

The approach that we have introduced in this paper can also be 
generalized to the more interesting case of proton Compton scattering 
PU| , where experimental data are availiable, while a direct comparison 
with perturbative QCD results based on leading power factorization 
theorems [^] is possible. Although the leading power result will, by 
definition, dominate for sufficiently large t, we expect the sum rule 
result to be a valuable tool in the interpretation of realistic experi- 
ments. Of particular interest will be the angular dependence of the 
sum rule result, which may distinguish between leading-power factor- 
ization theorems and the "Feynman mechanism" that dominates the 
sum rule spectral functions. This could help us to understand the 
roles of these two formalisms for describing hadronic elastic scattering 
at available energies. 

Acknowledgments. This work was supported in part by the Na- 
tional Science Foundation, grants PHY-9211367 and 91080541, by the 
U.S. Department of Energy under contract DE-AC05-84ER40150, and 
and by the Texas National Research Laboratory. 

A Analyticity and Infrared Safety 

To derive a sum rule for the Compton scattering amplitude, we must 
be able to prove, first, that the amplitude satisfies the dispersion re- 
lation, eq. (|I~9"D. In addition, we must check that for nonzero values of 
pj the amplitude is infrared safe; i.e. that its perturbative expansion 
is free of infrared divergences. Otherwise, the lowest-order term in 
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the expansion will be swamped by divergent higher-order corrections. 
The demonstrations of these two features of the fixed-angle Comp- 
ton amplitude are related, and we present them together. We will, of 
course, work in perturbation theory. 

Eq. (|19|) is consistent with singularities at positive values of p\ 
and/or p\, corresponding to normal and anomalous thresholds. Let us 
recall [|l(| that the condition for a singularity in a Feynman amplitude 
is that the momentum integrals be "pinched" between coalescing poles 
from the propagators of internal lines. If internal lines are not on-shell 
the integrand is finite; if the momentum contours are not pinched, we 
may deform the momentum integrals into another region where the 
integrands are finite. 

A useful way of characterizing points in momentum space where 
some set of lines go on-shell is through "reduced diagrams" , in which 
all lines that are off-shell are shrunk to points. Now suppose the 
resulting reduced diagram represents a point in momentum space at 
which the contours are truly pinched. Then a theorem due to Coleman 



and Norton p3| states that the diagram corresponds to a classical 
scattering process, in which all vertices correspond to points in space- 
time, connected by the motion of on-shell, free, classical particles, 
whose velocity is determined by their mass and momentum. This 
theorem turns out to be a very useful tool in the analysis that follows. 

Let us start by considering Compton scattering with space-like 
photons, qf < 0. We will relax this condition later. Now suppose p\ 
is either negative, or complex, but with an amplitude that obeys the 
restriction, 

b-|<A 2 , (55) 
with A 2 = (s + i)/4, as in eq. (15), so that 

Re{u) < 0. (56) 

We will be trying to show that for such a p\ the Compton amplitude 
has cut plane analyticity in the p\ plane when p\ also satisfies eq. (^5|) . 

Let us see what kind of physical-scattering reduced diagrams, which 
potentially give singularities, can include external currents with mo- 
menta in the ranges just described, all the time keeping s and t fixed 
and large. Because s is positive and t negative, we may assume that 
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the energies of p± ... q 2 remain positive. The first vertex in the re- 
duced diagram must therefore result from the scattering of p\ and q±, 
or of some subset of particles from one of their virtual states. When 
q\ < and p\ < or complex, however, neither the photon nor pi 
can decay into only on-shell physical particles as part of the physical 
process. Since t < 0, p\ also cannot decay by emitting the (possibly 
unphysical) momentum p 2 , and because u < (in the region we are 
considering), it cannot decay into a physical state by emitting q 2 , even 
if p\ is real. 

We may now take the limit q\ — > 0, which brings us up to, but not 
beyond, thresholds in the q\ channel. The virtual states into which the 
incoming photon can "decay" consist of collinear massless particles, 
which propagate together at the speed of light. In any space-time 
picture, they therefore act in the same manner as a single particle. 
The first vertex in the reduced diagram must therefore involve the 
direct scattering of q\ (or its collinear, lightlike decay products) and 
Pi at a point, and all the remaining particles must emerge from that 
point. Any situation such as this, in which all particles emerge from 
a single point in space-time, strongly reduces the number of possible 
reduced diagrams |24|. In fact, in this case, the only possible reduced 
diagrams corresponding to singular points involve no further "hard" 
scatterings with nonzero momentum transfer. All further interactions 
involve the local interactions of massive particles, relatively at rest, 
of exactly collinear-moving massless particles, and/or the emission 
and absorption of exactly zero-momentum particles. In our case, this 
means that we may only have thresholds when p\ = mi) 2 , for some 
set irti of the masses of virtual particles. Similarly, for q 2 — > we may 
have on-shell massless virtual particles moving exactly parallel to q 2 . 
We need not consider soft particles that connect lines moving parallel 
to p 2 with those moving parallel to q 2 , because, since the photon is 
color and charge neutral, its couplings to soft particles vanish on-shell 
Z5|. In addition, since we keep q 2 fixed, we only need to consider 



the singularities associated with p 2 , that is, its normal thresholds for 
p\ > 0. Clearly, analogous reasoning, working backwards from the 
final toward the initial state can be used to show that the amplitude 
has cut-plane analyticity in p\ for p\ fixed and off the real axis. This 
is the first result we set out to prove. 
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Given the preceeding arguments, the proof of infrared safety is not 
difficult. We are interested in the behavior of the amplitude around 
the contours at \pf\ = A 2 . In this region, we want to work in perturba- 
tion theory with massless internal lines. It is easy to see that the only 
reduced diagrams which correspond to points at which the momentum 
contours are pinched consist of a single hard scattering vertex, and of 
two "jet" subdiagrams, connected to the external photons, and each 
consisting entirely of lines collinear to that photon. There may also 
be zero momentum lines connecting to two jet subdiagrams and/or 
the hard vertex. Note that because we work to lowest order in QED, 
the jet subdiagrams are not photon self-energies, but must attach to 
the hard scattering by at least two finite-momentum lines each. The 
second result claimed above, the infrared finiteness of such an ampli- 
tude, then follows directly from the arguments given in ]2S| for e + e~ 
annihilation into photons. 



B Kinematics 

In this appendix we discuss the kinematics of the pion Compton scat- 
tering process in Fig. 3. 

Let qi and q 2 be the momenta of the incoming and outgoing pho- 
ton, respectively, which are assumed to be on-shell (g 2 = gf = 0). Let 
also pi and p 2 be the momenta of the incoming and outgoing pion. 
The external pion states are off-shell and are characterized by the 
invariants S\ = p\ and s 2 = p\- 

We define as usual the Mandelstam invariants 

s = (pi + <?i) 2 , t = (q 2 - qi) 2 , u = (p2 - <?i) 2 , (57) 
with the usual relation 

s + t + u = s 1 + s 2 . (58) 

We consider both s and t to be very large and in the physical region 
s > 0,t < 0. 

It is also convenient to introduce light cone variables for the pion 
and photon momenta as follows: 

Pi = pfv + PiV, p 2 =p%v + p2V 1 (59) 
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qi = qt v + Qi v + Qi-L, (60) 
v = ^(1,1, ± ), v = -j={1,-1,0 ± ), 

q 1± ■ n± = 0. (61) 

In the definitions of v and v, the order of components is (0,3,2,1). 
The momenta of the pions are, in this frame, purely longitudinal. We 
can choose a particular brick- wall frame in which pt = p 2 = Q- Then 
we have 

A covariant expression for the momentum scale Q can be found 
from the relation 



t=(P2-pif = s 1 + s 2 -2Q 2 -^. (63) 



We easily find, by solving for Q 2 , 



Q 2 = ^( Sl + s 2 -t + 5) = ^(s + u + 5) , (64) 



where we have defined 

5 = yJ( Sl + S2 - t) 2 - As lS2 = 2Q 2 - S 0. (65) 

In the scattering process at high energy and fixed angle, Q 2 is a large 
parameter. Note that in this notation t = —1Q 2 only when s\ = s 2 = 
0. In this frame we also easily get 

u = {p 2 -q l f = 2Q 2 -s+-^0. (66) 
Covariant expressions for qf and q 2 can also be found in the form 

qt = (67) 

22 



The magnitude of the transverse component q±± of the momentum of 
the incoming photon is then given by 



;_ J )( s _2Q2)(2Q2 s - SlS2 ) 

1,11 s 91 = 75^ ' (69) 

Similar formulas hold for the momentum of the outgoing photon: 

(2Q 2 -s 2 )(2Q 2 s-s lS2 ) 



4Q 3 S 
(s-2Q 2 )(2Q 2 - Sl ) 



C Computation of the Spectral Weights 



In this appendix, we discuss the evaluation of the integrals in eq. 
and their asymptotic behavior at high energy. The results for the 
crossed diagram are carried out analogously. 
We first note that we may use the identity 

2q 1 -k = s-s 1 -(p 1 -k + q l ) 2 (71) 



to rewrite eq. ( f47|) in terms of integrals for the fully on-shell box 
and triangle diagrams that are no more than quadratic in in the 
numerator. Thus, we have 

A (12W = 16(gi . pi + Ql . pa ) J d ± k J(4a) pi ) _ („ . fc )2 J 

+ 8 | d 4 k [1^ -( S - Sl )I^} 

x ((n ■ k)(n ■ pi + n ■ p 2 ) - 2(n • A;) 2 ) , (72) 

where 1^ and J^ 4a - ) are the integrands for the double discontinuities 
of the triangle and box diagrams, respectively, 

If) = 5 + (k 2 )5 + (( Pl - kf)5 + ((p 2 - kf) , 

r(4.) Mfc 2 M+((Pl-fc) 2 M+((P2-fc) 2 ) ^ 

s " (Pi-^ + gi) 2 ' 1 J 
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The relevant integrals are now 

= J lf ] k^k v (74) 

for the triangle and 

< = J ^ a) k" 
W£ u = J I^ a) k"k u (75) 

for the box. 

Before discussing the foregoing integrals, it is useful to describe 
the corresponding scalar integrals for the scalar triangle and the scalar 
box. For the triangle we have 

A 3 = J d 4 kS + (k 2 )S + (( Pl - k) 2 )5 + ((p 2 - kf) = ^, (76) 

where 5 has been defined in eq. (]65f) . The integral for A3 can be simply 
evaluated by going to the brick-wall frame defined above; 25 appears 
as the Jacobian of the transformation from the components fc M to the 
variables k 2 ,p 1 ■ k, and p 2 ■ k, leaving one trivial angular integral, 

with the angle in transverse space. Interestingly, the discontinuity 
on the cut is given by the same result for massive as for massless 
quarks. 

For the scalar box diagrams in Figs. 5(a) and 5(b), three internal 
lines are on shell, while the upper internal line is off-shell, 

A a {s 1 ,s 2 ,s,t) = d k — . 78 

J (p 1 -k + qi) z 

Three angular integrals are trivially done using the change of variables 
mentioned above, but the remaining angular integral is now nontrivial, 

A s a c ( Sl ,s 2 ,s,t) = j- r dJ \ -), (79) 

Ad Jo \s - si - 2(q 1L ■ k L - q L k±cos(p) J 



24 



where we define 

q 1L ■ k L = q+k~ + gf/c + , (80) 

which is a number that is fixed in the brick-wall frame according to 
the values given in Appendix B and the mass shell delta functions (see 
below). The angular integral is a standard one, 



o \A + Boost) ~ (A 2 - B*)V* • ( ] 

For the particular constants 

A = s - si - 2q 1L ■ k L = —j^—[s(s + u) - 2sis 2 ] , 

B = 2q ± k ± = 2 -^[s{-u) + s l s 2 ] 1 / 2 [s 1 s 2 ] 1 l\ (82) 



we find the following, rather unobvious, identity, 
1 1 



(83) 



(^2 - 5 2 )V2 ( [s _ Sl _ 2 q 1L ■ k L f - [2q ± k ± ] 2 y/ 2 s(-t) ' 

Thus, at lowest order, the box and crossed box are given by the re- 
markably simple forms, 

A - = -h • (84) 

where we have used the crossing relation between the two. 

In deriving the spectral weight for the box diagram with fermions, 
it is useful to know expressions for the components of in the brick- 
wall frame kinematics of Appendix B, 

_ s 2 (Q - sjpQ) 
k ~ 6 ' 

k- = ( 86 ) 
o 

as fixed by the delta functions in the integral. We notice that both 
components are small and 0(1/Q). Similarly, we have the simple 
expression 

k ± = — j 2 . (87) 
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Note that all the components of vanish in the limit Q — > oo, as is 
appropriate for the Feynman mechanism (Section 5). 

The integrals that define the spectral functions, eq. (fj7|), are eas- 
ily reduced in the brick- wall frame to the factor 1/45 times angular 



integrals which generalize eq. (^T 

/■2w 



C0S - ' (88) 



o \A + Bcos 

with a = 0, 1, 2 and with A and B given as above. For completeness, 
we list the elementary results, 

2tt 2tt<5 



2nA 2tt 4ir5 \f s(-u)sis 2 + s?s| 



(A 2 - B 2 ) 1 / 2 s(-t) 

h = 
h = 



B(A 2 -B 2 y/ 2 B s 2 (-t) 4Q 2 -2 Sl s 2 /s 
2nA 2 2nA 2tt5 s(s + w) - 2sis 2 



B 2 (A 2 - B 2 ) 1 / 2 B 2 s 2 (-t) 4Q 2 -2 Sl s 2 /s 



(89) 



We note that in the frame we work, A ^> B, and that while Jo and 
I 2 behave as 1/A for large A, Ii goes as 1/A 2 . In terms of these 
expressions, we find, 

7T 

= -^[s 2 (p 1 -p 2 - si)p^ + s x (pi -p 2 - S 2 )P2 ] , (90) 
= — (v»v u - Id^v 2 ) , (91) 

7T V 2 / 

where d 11 = d 22 = 1, with all other components zero, and 

^=^ + ^S&*- (92) 



7T 

' SJ _ -(12^2^2 + (/0 - I 2 )V^l) • (93) 



45 3 
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The leading high-energy behavior of /\(^ 2 )p ert in eq. ( f4T|) comes 
entirely from the term 

A (12)pert = 8 ( s _ u )( n . u , a )( n .p 1 ) + ... 

~ 4 : (s-u)(^^(s 1 + s 2 yj(n-pi)(n-p 2 ) , (94) 



where we have used (n ■ pi) = (n ■ p 2 ). Note that by eq. (|5C)D , the 
factor (n ■ pi) (n ■ p 2 ) cancels in the sum rule. Similarly, we find for the 
crossed diagram Aj, 12)per *: 

Af )pert ~ 4(« - s) + (n • px)(n • p 2 ) . (95) 

Finally, using eqs. ([46]) and (0), we get the asymptotic spectral func- 
tion quoted in eq. (|5T|) . 
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Figure Captions 



1. The integration contours for the Borel transform in the pf plane, 
eq. ©. 

2. Lowest order perturbative contribution to the pion form factor. 

3. The pion Compton scattering amplitude. 

4. The integration contours for the Compton scattering amplitude. 
Here 7 is the dispersion contour of eq. (||), and T the transform contour 
ofeq. 

5. Lowest order contributions to the spectral weight for Compton 
scattering. 
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